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Abstract 

We study the pairing instability and mechanical collapse of a dilute homo- 
geneous bose gas with an attractive interaction. The pairing phase is found 
to be a saddle point and unstable against pairing fluctuations. This pairing 
saddle point exists above a critical temperature. Below this critical tempera- 
ture, the system is totally unstable in the pairing channel. Thus the system 
could collapse in the pairing channel in addition to mechanical collapse. The 
critical temperatures of pairing instability and mechanical collapse are higher 
than the BEC temperature of an ideal bose gas with the same density. When 
fluctuations are taken into account, we find that the critical temperature of 
mechanical collapse is even higher. The difference between the collapse tem- 
perature and the BEC temperature is proportional to (nla^P)^''^, where n is 
the density and is the scattering length. 
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The system of Bose gas with attractive interaction has received a lot of attention recently 
due to the success in cooling ^Li [1] and ^^Rb [2] systems. In these trapped systems, Bose 
Einstein condensation can be realized at zero temperature if the number of atoms in the 
condensate is below a critical value [3], because of a finite size effect. In a ^Li system, 
the condensate experiences repetitive growth and collapse as the system is cooled down [1] . 
In a ^^Rb system, the interaction can be varied from being repulsive to being attractive 
by tuning a Feshbach resonance. During this process, the system loses particles until the 
critical number is reached [2]. 

A homogeneous Bose gas with attractive interaction is generally unstable at zero tem- 
perature because the system can always lower energy by increasing density. At high tem- 
peratures entropy stabilizes the gaseous phase. There are several scenarios as to how the 
system becomes unstable when the temperature is lowered. The simplest picture is that the 
system collapses within the normal gaseous phase, which is supported by several theoretical 
studies [5] . This collapse is mechanical collapse, where the compressibihty becomes negative 
and the system can lower energy by separating into a low density phase and a high density 
phase. Another possibility is that the system goes into the Bose-condensed phase and then 
collapses in the condensed phase. However the condensation transition temperature is found 
to be lower than that of mechanical collapse [5]. In addition to these scenarios, Evans and 
Imry proposed that there is a pairing phase similar to the BCS phase in superconductors at 
low temperature. However it was also concluded in the previous theoretical studies [5] that 
this pairing phase transition temperature is lower than the mechanical collapse temperature 
of the normal phase. 

It was not fully addressed in the previous studies whether the pairing phase could exist 
as a metastable phase or not. In this paper, we examine its stability in detail. In addition, 
we go beyond the Hartree-Fock approximation to see how the temperature of mechanical 
collapse is affected by fluctuations. 

The interaction between atoms is singular at short distances. However, as for other cases 
in condensed matter and elementary particle physics, at low energies the dilute Bose gas 
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system can be well described by the Hamiltonian with an effective contact interaction [6] with 
the coupling constant proportional to the scattering length a^, up to a momentum transfer 
which is determined by the effective range of the interaction. The unphysical ultraviolet 
divergences in this approach can be subtracted by the pseudopotential method [6]. This 
renormalization scheme is justified in the case of a dilute Bose gas with weak interaction at 
zero temperature, where it can be shown that the divergences come from double counting 
of certain diagrams and the dominant perturbations come from the repetitive scatterings of 
two particles. 

At finite temperature when there is no condensate, the medium effect becomes important 
and the coupling constant could depend on other quantities in addition to Og. To avoid 
this complexity, here we use a simple model with a contact interaction which has coupling 
constant g and an explicit cutoff A in the wave- vector space. The Hamiltonian of the system 
is given by 

H = ^VV'^ • Vip + #^ V- (1) 

A length scale a can be defined from the coupling constant, a = m/ {4:TTh'^g), but a is not the 
scattering length Ug. By solving the two-body scattering problem, we obtain the relation 
between a and ag, = a/{l + 2aA/7r). The cutoff A is proportional to l/^\as\rs, where 
Ts is the effective range. Here we consider the case of a weakly attractive interaction with 
a < 0, — 2aA/7r < 1, and \na^\ <S 1, where n is the density. 

To study this system, we use the Peierls variational method [7] which gives an upper 
bound for the free energy or grand thermodynamic potential of a quantum system. To discuss 
both the normal phase and the pairing phase, we assume the quasiparticles are superpositions 
of particles and holes, and they have certain occupation numbers at finite temperature. In 
momentum space, the density matrix is given by < V'kV'k >= (/fe + 1/2) cosh^^ — 1/2, and 
< V'kV'-k >= {fk + l/2)sinh^fe, where cosh^^ and sinh^fe are coherence factors and are 
chosen to be real. The occupation number is given by fk- 

Within this approach, the expectation value of the grand potential is given by 

3 



^(m, {fk}) = J2i^k - [{h + 1/2) cosh - 1/2] + + 1/2) sinh O^f 

+|^{E[(A + 1/2) cosh^fe - 1/2] }2 - T5][(/fe + 1) ln(/fe + 1) - /feln/^]. 
^ k k 

The normal phase and the pairing phase are given by the saddle points of the variational 

space, dQ/dfk — 0, dQ/d9k — 0. The saddle point equations can be further written as 

tanh^fe = , (2) 

efe - A - /i 

A = ^coth§-l (3) 



where 



^= -^E(/fc + 1/2) sinh ^fc, 
^ k 

^=-T7E[(/fc + V2)cosh^fe-l/2], 

k 



Ek — ^Jfil — A^, and rjk — — A — ji. The pairing phase corresponds to the solution with 
finite A and the normal phase corresponds to the solution with A = 0. The transition from 
normal phase to pairing phase occurs when A approaches zero. 

To simplify the discussion without compromising the saddle point structure, we look at a 
more restricted variational space defined by Eq.(2) and Eq.(3) with A and A as variational 
parameters. In terms of these new variational parameters, the grand potential is given by 



(1) 



where the parameters A, A are defined by 



A(A, A) = y (/fe + -) sinh = - A y A coth ^, 



and 



A(A, A) = -2^ y [(A + -) cosh 6^--] = -- y [— coth ^ - 1]. 
The new saddle point equations are given by A = A, A = A. 
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The stable phases correspond to the local minimum points of the grand potential. If the 
pairing phase is stable, one of the necessary conditions is that it should be stable against 
fluctuations in pairing, which is given by 



- . >0. (5) 



9A2 

Since the transition into the pairing phase would be a second-order phase transition, at the 
critical temperature Tp, 

and the stabihty condition Eq.(5) becomes 



9A4 

However, we find that the opposite inequality holds 



> 0. (7) 
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-Ett-3(1 + 2/.^)<0. (8) 



A=A, A=o 2g ^ 27]l 

It clearly shows that the pairing phase is generally unstable. The instability of the pairing 
phase can also be illustrated by plotting the profile of the grand potential as shown in 
Fig.(l). The pairing phase is located at a saddle point and the normal phase is located at 
a minimum point. The pairing saddle point is a maximum in the direction of the normal 
phase and a minimum in the perpendicular direction. 

This instability can be studied in more detail by constructing an effective Lagrangian 
near the critical temperature Tp based on Landau's theory of phase transitions 

£(A) = 0( A, A) - 0(0, A) f« aA^ + /?A^ (9) 

where a = ^U=a,a=o and /3 = aSl^ Ia=a,a=o> < from Eq.(8). When T > Tp, 
a > 0, there is a pairing saddle-point solution A^ = — q;/(2/3). When T < Tp, a < 0, the 
normal phase becomes unstable in the pairing channel and there is no pairing saddle-point 
solution. It is clear that the system is unstable in the pairing channel when the temperature 



is below Tp. Even above Tp, the normal phase can become unstable in the pairing channel 
by tunneling into the region beyond the pairing saddle point if any place in that region has 
lower energy than the normal phase. The pairing saddle point serves as an energy barrier 
in this case. However, in our limited numerical calculations, we have not been able to 
determine the condition for the normal phase to become metastable. When the temperature 
is far above Tp, the pairing saddle point disappears. In this case the effective Lagrangian 
description breaks down and a better method is needed to find the general expression of the 
critical temperature below which the pairing saddle point appears. 

The critical temperature of pairing instability Tp is given by the following equation 

l^_^/^J_coth^, (10) 

where r]k = + 2gn — /j,. This equation has solution only when the temperature is close to 
the ideal BEC temperature Tq. In this case it can be further simplified. 



1 + 

and 



2aA ^ r d'k Tp 



2aA Tpa 2m 
1 + «-^/_. 12 



From the density equation, we obtain 



where 



""-iL^^fl'^^-'*^- '''' 



_ 2TTh , n , 2 

Tn= -T^ 3 



^2' 

and C,{x) is the Riemann-Zeta function. Therefore the solution of Eq.(lO) is given by 

Tp 27r 1 1 I Stt 1 , , 

To 3Ct(i)(l + ^) 3Ct(i) 

In addition to the pairing instability, mechanical collapse can also occur in the normal 

phase as discussed in several references [5]. It is important to find out which instability 
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takes place first. The mechanical collapse happens when the compressibility goes to zero 
d/i/dn — 0, which can be more explicitly written as 

/7^csch^^ = l, (15) 



2TcJ (27r)3 2Tc 

where Tc is the temperature of mechanical collapse. Similar to Eq.(lO), Eq.(15) can only 
be satisfied when Tc is close to Tq, and can be simplified as 

l--af^^ (16) 

and 



The solution is given by 



(18) 

To 3Ct(|) ' ' ^ ^ 

Similar results have been obtained in the previous studies [5]. 

By comparing the two critical temperatures Eq.(14) and Eq.(18), we find that when 
-AaA/n > 1, Tp > Tc; when -AaA/n < 1, Tp < Tc- However, this result is obtained 
within the variational approach, which is equivalent to the Hartree-Fock approximation. To 
obtain a better quantitative result, fluctuations need to be considered. 

There are infrared divergences in all orders of perturbation at Tq. To obtain the correct 
renormalization near Tq, summing the leading divergent terms at each order is necessary, 
which can be effectively performed in a self-consistent approach. Here we consider the 
simplest one-loop renormalization with self consistency. In the Hartree-Fock approximation, 
the self-energy in the normal phase is simply given by 2gn. When fluctuations are taken 
into account, the coupling constant g is renormalized. In one-loop order, g is re normalized 
by particle-particle scattering. It is well known that at zero temperature, such scatterings 
renormalize g to Anfi^as/m. In this approach, the renormalized coupling constant is given 
by the i-matrix 



t{q, n)=g-Tj2 77^s9t{ci, l^)G(k, uj)G{q -k,n-uj), (19) 

where G(k, a;) = l/(icc; — S(k, cu) — + /x) is the finite-temperature Green's function. Under 
this approximation, the new self-energy is given by 

E(k,a;) = 2rX: / ^t(q, 0)G(k - q, a; - O). (20) 

The infrared behavior of the integrals in Eq.(19) and Eq.(20) is dominated by E(0, 0) — 
which is small and finite at the collapse temperature. The A;-dependence and a;-dependence 
of E just provides a renormalization factor to the propagator in the leading order. So here 
we only consider only the constant part of self-energy, 

E = E(0,0) ^ 2t(0,0)n, 

and get 

where r;^ = -|- E — /x. Using the same technique in solving Eq.(lO), we obtain t(0,0) in 



terms of E 



i(0, 0) « — ^- — ; . (22) 



Therefore the self-energy E is given by the self-consistent equation 



-1 I 2aA I Ta I 2m 



(23) 



The condition of mechanical collapse dji/dn — Q now leads to the following equation 

Tm (9E / 2m ,„ , 

(25) 



where from Eq.(23) we obtain 
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The solution of Eq.(25) is approximately given by 



(26) 



(1+ (27) 



^2m(S - 1^) 
Close to To, the density equation yields 

where 5T — T—Tq. Combining Eq.(27) with Eq.(28), we obtain the new critical temperature 
of mechanical collapse 

T'c ,^,(27r|a>t; 



In 2 
3 



(29) 



^0 3Ci(|) 

The difference between Tq and Tq is now proportional to (n|as|^)^/^. In the weakly- 
interacting limit, it is much bigger than n^/^la^l which is proportional to the temperature 
difference in Hartree-Fock approximation. The strong renormalization of the collapse tem- 
perature is a result of the enhancement of the Hartree-Fock attractive interaction due to 
fluctuations. In contrast, the critical temperature of pairing instability is not renormalized 
by the i-matrix. We have used the one-loop and constant self-energy approximations which 
may limit the accuracy of the numerical coefficient. Future works with better approximations 
may improve the numerical prefactor and provide an estimate of the next order term. 

Our result shows that the collapse temperature is higher than the temperature of pair- 
ing instability. However the pairing saddle point exists above the temperature of pairing 
instability. The saddle point may still be important to collapse dynamics if it appears above 
the collapse temperature. In this case, if the normal phase is only a local minimum and 
not a global minimum of the grand potential, there is a finite probability for the system to 
pass the pairing saddle point and tunnel into the unstable region. Discussing the dynami- 
cal process of collapse is beyond the scope of this paper. However, within this variational 
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framework, the collapse process can probably be studied by the general method of analyzing 
the stabihty of metastable states [8] in the future. 

We would like to thank the Aspen Center for Physics where some parts of this work 
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No. DMR 98-15932. 
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